Quantum critical fluctuations in disordered rf-wave superconductors 
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Quasiparticles in the cuprates appear to be subject to anomalously strong inelastic damping 
mechanisms. To explain the phenomenon, Sachdev and collaborators recently proposed to couple 
the system to a critically fluctuating order parameter mode of either id^y- or is-symmetry. Motivated 
by the observation that the energies relevant for the dynamics of this mode are comparable to the 
scattering rate induced by even moderate impurity concentrations, we here generalize the approach 
to the presence of static disorder. In the id-case, we find that the coupling to disorder renders 
the order parameter dynamics diffusive but otherwise leaves much of the phenomenology observed 
in the clean case intact. In contrast, the interplay of impurity scattering and order parameter 
fluctuations of is-symmetry entails the formation of a secondary superconductor transition, with a 
critical temperature exponentially sensitive to the disorder concentration. 
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Recent experimental work shows that essential aspects 
of the low-temperature behaviour of cuprate supercon- 
ductors lie beyond the scope of orthodox BCS-type for- 
mulations in terms of a non-interacting quasiparticle 
(QP) system. Notably, high resolution angle-resolved 
photoemission (ARPES) experiments indicate a QP 
relaxation rate k ~ jej proportional to the excitation en- 
ergy (for temperatures T < |e|), while the standard BCS 
picture predicts k, ~ \ef Q|. The much weaker ^ |e 
scaling most likely represents a correlation effect. 

Arguing that a qualitative modification of the |ep-law 
required coupling of the QP system to some kind of col- 
lective mode, Sachdev and collaborators 1^ recently pro- 
posed a dynamical generalization of the BCS order pa- 
rameter field, viz. 



A(x,i) = Ao(x) +iAi(x,i), 



(1) 



where Aq is the background field of dxi-yi angular 
momentum dependence and iAi = iA^^'^^y a time- 
dependent component of either s- or d^jy-wave symmetry. 

The basic idea behind the proposal (|^) is that some- 
where in close vicinity to the phase diagram of the sys- 
tem there might be a quantum critical point towards 
formation of a stable iAi order parameter amplitude (a 
d + is/dxy state.) Although the actual transition point 
may lie outside the accessible parameter space, its prox- 
imity will manifest itself through dynamical order param- 
eter fluctuations, as described by (|l]). 

The coupling of the collective mode Ai{t) to the rela- 
tivistic QPs strongly alters the low-energy phenomenol- 
ogy of the system ||, Q . Specifically, (i) scattering off the 
time-dependent order parameter fluctuations enhances 
the QP relaxation to a linear rate k ~ |e|. This amplifica- 
tion mechanism appears to be insensitive to microscopic 
details such as the ultraviolet structure of the QP disper- 
sion relation. In particular, (ii) no qualitative differences 



between s- and (ia;j,-scattering exist. The coupling of the 
QPs to Ai leads (iii) to the formation of long-ranged 
collective fluctuations which should |^ be observable by 
neutron or Raman scattering. However, (iv) if sufficiently 
weak, the attractive interaction mediated by Ai does not 
open a QP gap, i.e., it does not drive the system into a 
secondary superconductor instability. 

It is the purpose of this Letter to explore whether ele- 
ments of the phenomenology outlined above survive gen- 
eralization to the presence of static disorder. In view 
of the fact that all features (i)-(iv) rely on a conspiracy 
of order parameter fluctuations and the Dirac spectrum 
of the QPs - the latter being thoroughly corrupted by 
impurity scattering - one might expect the answer to 
be largely negative. However, as we are going to show 
below this presumption is overly conservative: owing to 
a decoupling of the order parameter from the impurity 
scattering vertex (ic?-case) and/or the renormalization of 
the polarization operator by infrared-singular diffusion 
modes (is-case), the disordered system, too, turns out to 
be strongly affected by fluctuations of Ai{t). 

To summarize our main findings, (i') notwithstand- 
ing the presence of disorder, the coupling of QPs to a 
dynamical order parameter of id-symmetry continues to 
induce an inelastic relaxation rate nin However, 
this contribution sits now superimposed on a structure- 
less elastic background KcI ~ t~^, where r denotes the 
impurity scattering time. Further, at low frequencies, 
e| < T^^, a logarithmically singular third contribution 
Ksing ^ lii(T|e|) to the relaxation rate k = Kin + + Ksing 
signals the formation of long-ranged diffusive excitations 
which eventually (for asymptotically low excitation ener- 
gies/temperatures) fully absorb the spectral weight car- 
ried by individual QP states, (ii') Unlike in the clean 
case, the system is no longer indifferent to the symmetry 
of Ai(t). While in the id-case, QPs remain (marginally) 
stable objects and (iii') the order parameter continues 
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to exhibit fluctuation behaviour similar to that found in 
Refs. H), a conspiracy of impurity scattering and dy- 
namical fluctuations of is-symmetry (iv') drives the sys- 
tem into a secondary superconductor transition. More 
specifically, at the critical temperature 



(2) 



the system enters a d + is superconducting state and the 
QPs become gapped. Here loq = min[ti;£), 1/r], where 
is a phenomenological parameter setting the maxi- 
mum (Debye) frequency of the order parameter fluctu- 
ations, 1/f denotes the amplitude of the order parame- 
ter fluctuations, and i^q the residual low-energy density 
of states (DoS) induced by impurity scattering. Notice 
that through the dependence of the DoS on thc impu- 
rity concentration, vo 1/r, the transition temperature 
varies exponentially with the disorder strength. Since 
the latter can be deliberately changed by doping, the ob- 
servability of a disorder-dependent BCS superconductor 
transition at low temperatures appears to impose a rel- 
atively straightforward criterion testing the presence of 
dynamical is-order parameter fluctuations in the system. 

To prepare our discussion of the interplay order param- 
eter fluctuations/disorder, let us briefly recapitulate the 
physics of the clean system. Denoting the two-component 
QP field by -(/'"i where a = 1,1,2,2 enumerates the 
low-energy nodes in the Brillouin zone, and the time- 
dependent amplitude of the order parameter by (t>{x), the 
low-energy physics of the system is described by a com- 
posite action S[(j), ip] = S[ip] + S[(j)] + Sc[(j>, tp]- Here, x — 
(t, x)"^ is a three-component space-time argument and 
Slip] = J d^x 1/;°^ [dr + tof(s?7" Vi^i + s^a2d2)] ip" the 
Dirac action of the unperturbed QPs where v-p is the 
Fermi velocity, 7 = i/|Ao|, t the tight binding energy, 
s° are sign factors depending on the node index, and (jj 
Pauli matrices in particle/hole space (5|. Further, 



S[cj>] 
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4!' 



(3) 

is a 0^-type action controlling the fluctuation behaviour 
of the order parameter field, (c, r, and u are phenomeno- 
logical constants, where r = marks the position of the 
quantum critical point into the dx2_y2 + is/dxy phase.) 
Finally, the coupling between the QP fields and the col- 
lective mode (j) is mediated by the scattering vertex 



V'] = A / d^x V'^^FVsV" 4>, 



(4) 



where A is a coupling constant, and F° a symmetry 
factor discriminating between the two different types 
of order parameters: F" = 1 (F" = (— )°) for the 
case of is {id) symmetry. Perhaps thc most direct 
way to understand the consequences of this coupling 
mechanism is to integrate over the QP fields. This 



leads to an effective (/)-action S'cffi^] = S' + <^S' 
where 5S[4>] = — ln(exp(— S'c[</), ■0])) measures the ten- 
dency of the system to create/annihilate a Cooper pair 
in response to fiuctuations of the order parameter and 
(...) = / P(?/'^V)exp(-S'[V'^ ?/'])(. ..) is the functional 
average over the QP action. For sufficiently weak cou- 
pling, (55[(/)] can be approximated by its second order ex- 
pansion in (j) (an RPA type approximation), i.e., 5S[<f>\ — 
^ J d^xd^x' (j}{x)xix,x')(j}{x'), where 

X{x,x') = tr [FV3e'^'^'(:r,x')F"'(T3e'^'"(x',x)] (5) 

is thc Cooper pair propagator, 'tr' stands for a trace over 
particle/hole indices, and fj"" {x,x') — {ip"-{x)ip'''^ (x')) is 
the Gorkov Green function. 

Evaluating x for the clean case (in which Q becomes 
diagonal in momentum and nodal space), one finds 

SS[^] ^ - [d\ \cl>k\^ \Ci - C2{u;l + k2)i/2l ^ (6) 



where J d^k = Tj^uj j denotes a summation over 
the three-momentum k = {ujm,'^)'^ , and ujm is a bosonic 
Matsubara frequency. While the constant contribution 
Ci can be absorbed into a shifted value of the param- 
eter r, the non-analytic operator ~ (k^ -I- a;^)^/^^^ - 
arising from the infrared singularity of the Dirac opera- 
tor ~ dominantly influences the scaling behaviour of the 
model, and holds responsible for the anomalously strong 
QP relaxation rate. 

As a corollary we remark that (ii') the structures above 
are insensitive to the symmetry of Ai(t). The differ- 
ence between an s- and a d-scattering amplitude, re- 
spectively, is that the former is constant while the lat- 
ter changes sign from node to node. However, being 
quadratic in the scattering vertex, the pair susceptibility 
operator (and therefore the induced action 5S) is oblivi- 
ous to this difference. Further, (iv') the rough estimate, 
Ci A^ dev{e)/e implies that the same linearity of 
the DoS, i/(e) ~ |e|, that led to the non-analytic struc- 
ture of 5S excludes the formation of a Cooper instability 
in the 0-action. In contrast, coupling to a fictitious 
QP species with an extended Fermi surface, v{e) = I'o, 
would lead to Ci(T) ~ X^vo \n{ujQ/T), i.e., the formation 
of a superconductor transition at a critical temperature 
Tc ~ Wo exp(— f/i'o)j where f = r/A^. 

How then do these structures change in the presence of 
disorder? Before addressing this question in quantitative 
detail, let us briefly summarize a number of key elements 
entering the theory of the weakly disordered |^ system: 
first, the picture developed in Refs. ||, || essentially re- 
lies on the linear singularity of the low-energy DoS. In 
contrast, disorder scattering leads to a randomization of 
spectral structures over scales ~ and, therefore, to a 
smearing of the low-energy cusp of the Dirac spectrum. 
Relatedly, the QP dynamics crosses over from ballistic 



at |e| > to diffusive at |e| < t~^. As relativistic QP 
dynamics is an essential input to the theory ^ one 
may expect |e| ^ r^^ to be a lower bound for its appli- 
cability. Finally, and in a way to be discussed in more 
detail below, the relaxation time (as well as DoS, spin- 
and thermal conductance) are affected by mechanisms 
of quantum interference, similar to the weak localization 
corrections known from the physics of normal metals. 

To quantitatively explore the ramification of these 
mechanisms in the present context ^ ^ , we consider 
the configurational average of the susceptibility operator 
(x(x, x'))dis- The most basic effect of impurity scattering 
will be that the two Green functions appearing in ac- 
quire a self energy, C/"^(ie„) t/"^(je„C) + iKci sgn(e„), 
where the real constants Kc\ — (2r)^^ and C are obtained 
for |e| < from (depending on the microscopic real- 
ization of the disorder) self-consistent Born or T-matrix 
perturbation theory 

We next observe that in the is-case the inclusion of 
a Green function self energy is partially compensated 
for by vertex corrections. Evaluating the quasiparti- 
cle/order parameter vertex within a ladder approxima- 
tion shown in Fig. |^ (the latter stabilized by the pa- 
rameter 7 1), one finds that the vertex is coupled 
to a mode X'„„'(q) = (-Dq^ -I- |e„| + |e„'|) similar to 
the singular 'diffuson' and 'Cooperon' modes of disor- 
dered metals. Here, D — {j + 7~"^)/(7r^i/o) is the diffu- 
sion constant of the d-wave superconductor and e„,„' are 
fermionic Matsubara frequencies. Including this mode in 
the calculation of the polarization operator, we find that 



(7) 



where, apart from a T-independent shift Sr {< r), the 
constant term Ci(T) = A^i'o ln(a;o/T') + Sr now displays 
the characteristic low-temperature profile of s-wave su- 
perconductors, and the ellipses stand for operators in- 
volving derivatives. The low-T singularity of Ci implies 
that at the critical temperature (^ the (/(-action becomes 
unstable, and a BCS transition into a phase with finite 
expectation value of the order parameter ((0) ^ 0) and 
gapped QP states occurs. The only significant difference 
to the BGS transition in conventional disordered metals 
is that presently the low-energy DoS - an essential 
ingredient to the buildup of a BCS instability - is gener- 
ated by impurity scattering (rather than by the existence 
of an extended Fermi surface). I.e., unlike in normal met- 
als, the critical temperature is exponentially sensitive to 
the disorder concentration. 

The key phenomenon prescribing much of the physics 
of the complementary id-case is that here the suscep- 
tibility operator is not affected by vertex corrections. 
Formally, this follows from the fact that the mode V 
obtained by summing an impurity ladder (cf. Fig. ^ is 
isotropic in node space. As a consequence, the two nodal 
summations over the symmetry factors F" = (— )" at the 
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FIG. 1: Top: structure of the polarization operator in the 
clean, disordered id, and disordered is case, respectively. Bot- 
tom: effective order parameter propagator. 

vertices decouple and annihilate all contributions to the 
susceptibility operator, safe for the 'bare bubble' without 
vertex corrections. Heuristically, this mechanism can be 
understood by noting that multiple impurity scattering 
occurring in transit between two consecutive interactions 
with the order parameter field distributes the QP ampli- 
tude homogeneously over all nodes. The vanishing of the 
nodal average of A'*^« then implies the absence of multi- 
ple scattering contributions to the pair propagator. 

Evaluating Eq. (j|) for the self-energy decorated Green 
functions, we find 



J d^k\cl>k\^ [C, ~ C2 {\uJrn\ + Dk^ 



(8) 



with the (now again finite at T 0) constant Ci ^ 
X'^uJofi^j^DT'C), where f{x) decreases from 1 at a; 1 to 
X at X — > 0, and C2 A^. Importantly, a non-analytic 
contribution ~ \ujm\(f>^ survives the blurring of the low- 
energy spectrum by the imaginary part of the self-energy. 
The continued presence of an operator of engineering di- 
mension 1 indicates that the collective mode displays 
fiuctuation behaviour similar to that in the clean case. 
Relatedly, a straightforward calculation along the lines 
of Ref. |4| shows that in the vicinity of the instability 
(r ^ C2 e|) Kin ~ as in the absence of disorder, while 
for r ^ C2|e| we obtain Fermi liquid behaviour Km e^- 

Having identified the two terms Kin and k^i, we proceed 
to discuss the contribution of disorder-induced quan- 
tum interference to the QP self energy. To leading or- 
der in an expansion in the small parameter g"^, where 
g = Di/Q, the soft modes V couple to the self-energy op- 
erator through the diagram shown in Fig. ^j. Physically 
(cf. Fig. H), this diagram describes the self-interference 
of a QP traversing a scattering path twice. In close anal- 
ogy to the weak localization corrections to the conduc- 
tance of normal metals, these processes lead to a negative 
logarithmic correction Ksing(e) ~ (7~-^ln(T|e|) to the self 
energy. In passing we note that the singularity of the in- 
terference correction to the DoS entails the existence of 
a crossover scale uj* ^ exp(— g), below which a pro- 
liferation of diffusion modes holds responsible for both, 
Anderson localization of the QP states and a linear van- 
ishing of their DoS, u[e) - vo\A/^* % HI- 

While the logarithmic singularity of the self energy 
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be made visible experimentally. 
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FIG. 2: Top: quantum interference correction to the self en- 
ergy (left) and corresponding Feynman path (right). Bottom: 
renormalization of the interference process by interactions. 

is primarily caused by disorder scattering, the presence 
of a dynamical id order parameter necessitates to re- 
place the bare diffusion mode by the 'interacting' mode 
shown in Fig. |^ (bottom). However, the quantitative 
summation over these self-interaction processes merely 
leads to an inessential renormalization of the prefactor, 
g^^ — > g~^{l + F), where F <g; 1 is an interaction- 
dependent constant, playing a role similar to the Fermi- 
liquid constants in normal metals. 

Similar logarithmic corrections, weakly renormalized 
by QP interactions, are found for the (spin-) conduc- 
tance and the DoS jl^, p^. These findings conform with 
the RG analysis of Ref. [|10| which showed (albeit for the 
case of globally broken time-reversal invariance - sym- 
metry class C in the notation of Ref. while the 
present problem falls into the more complex symmetry 
class CI) that interactions only weakly affect the low- 
energy disorder-generated interference phenomena . 

Summarizing, we have explored the interplay of dis- 
order scattering and critical fluctuations of an is- or id- 
order parameter component in two-dimensional cuprate 
superconductors. In the zs-case, the most prominent ef- 
fect caused by disorder is a secondary BCS transition 
whose critical temperature sensitively depends on the im- 
purity concentration. In the complementary id-case, the 
addition of disorder to the system leaves much of the 
phenomenology derived in Refs. [|[ Q intact. Notably, 
for temperatures T < |e|, the characteristic singularity 
Kin ~ |e| pertains to the disordered case. Heuristically, 
these phenomena can be understood by noting that in the 
presence of disorder (and for mass parameters fine tuned 
to criticality) the collective mode corresponding to the 
order parameter amplitude diffuses and, therefore, con- 
tinues to act as a strong agent of QP scattering. Finally, 
we included the (essentially disorder-generated) forma- 
tion of logarithmic singularities at low energies into our 
discussion of the relaxation rate. Given that even mod- 
erate impurity concentrations give rise to scattering 
times comparable to the energies (O(20K)) estimated ||] 
as characteristic for the fluctuations of Ai, we believe 
that, for sufficiently low temperatures, our findings may 



[1] 
[2] 

[3] 



[4] 
[5] 
[6] 

[7] 

[8] 

[9] 



[10] 

[11] 
[12] 



[13] 
[14] 



[15] 
[16] 

[17] 



Also at A.F. loffe Physico- Technical Institute, 194021 

St. Petersburg, Russia. 

T. Valla et al, Science 285, 2110 (1999). 

S.M. Quinlan, D.J. Scalapino, and N. Bulut, Phys. Rev. 

B 49, 1470 (1994). 

See, e.g., M. Vojta, Y. Zhang, and S. Sachdev, Phys. Rev. 
B 62, 6721 (2000); M. Vojta and S. Sachdev, Adv. Solid 
State Phys. 41, 329 (2001). 

D.V. Khveshchenko and J. Paaske, Phys. Rev. Lett. 86, 
4672 (2001). 

See A. Altland, B.D. Simons, and M.R. Zirnbauer, Phys. 
Rep. 359, 283 (2002) for review. 

We assume Ao to be so large that the pair breaking effect 
of disorder can be neglected. Furthermore, we do not 
consider intrinsic randomness in the parameter r Q as it 
does not qualitatively change our results. 
See, e.g., S. Sachdev, Quantum phase transitions (Cam- 
bridge University Press, Cambridge, 1999). 
S. Vishveshwara, T. Senthil, and M.P.A. Fisher, Phys. 
Rev. B 61, 6966 (2000); M. Fabrizio, L. Dell'Anna, and 
C. Castellani, Phys. Rev. Lett. 88, 76603 (2002). 
For a previous discussion of the effect of A* fluctua- 
tions on the disordered QPs, see D.V. Khveshchenko, 
A.G. Yashenkin, and I.V. Gornyi, Phys. Rev. Lett. 86, 
4668 (2001). However, due to an inconsistent treatment of 
the interaction propagator in this paper, the BCS insta- 
bility was overlooked and, instead. In ln(T|e|)-type quan- 
tum corrections to the DoS and the conductance were 
found for f/^'o < 1. 

M. Jeng, A.W.W. Ludwig, T. Senth il. and C. Chamon 



BuU. Am. Phys. Soc. 46, 231 (2001); |cond-mat/0 112044 . 
PA. Lee, Phys. Rev. Lett. 71, 1887 (1993); A.G. 
Yashenkin et al, Phys. Rev. Lett. 86, 5982 (2001). 
In the is-case, for to* ~ Tc, the BCS instability may inter- 
fere with the formation of strong localization phenomena. 
As oj* and Tc depend on fundamentally different input 
parameters, either case ui* < Tc ot lj* > Tc is conceiv- 
able |l|]. 

J.S. Meyer, I.V. Gornyi, and A. Altland, to be published. 
For a related study on the impact of critical fluctuations 
on the residual conductivity o f two-dimensional sy stems, 
see Y.B. Kim and A.J. Millis, |cond-mat/021026C| ! 
M.R. Zirnbauer, J. Math. Phys. 37, 4986 (1996); A. Ah- 
land and M.R. Zirnbauer, Phys. Rev. B 55, 1142 (1997). 
While these statements apply for a generic form of the 
interaction (r 7^ 0) , fine tuning the collective mode into a 
critical regime will lead to more pronounced interaction- 
generated phenomena [^ in both, is and id, cases. 
T.R. Chien, Z.Z. Wang, and N.P. Ong, Phys. Rev. Lett. 
67, 2088 (1991). 



